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It is shown that it is possible to bosonize fermions in any number of dimensions using the 
hydrodynamic variables, namely the velocity potential and density. The slow part of the Fermi 
field is defined irrespective of dimensionality and the commutators of this field with currents 
and densities are exponentiated using the velocity potential as conjugate to the density. An 
action in terms of these canonical bosonic variables is proposed that reproduces the correct 
current and density correlations. This formalism in one dimension is shown to be equivalent 
to the Tomonaga-Luttinger approach as it leads to the same propagator and exponents. We 
compute the one-particle properties of a spinless homogeneous Fermi system in two spatial 
dimensions with long-range gauge interactions and highlight the metal-insulator transition in the 
system. A general formula for the generating function of density correlations is derived that is 
valid beyond the random phase approximation. Finally, we write down a formula for the annihi- 
lation operator in momentum space directly in terms of number conserving products of Fermi fields. 

PACS numbers: 71.10.Pm,73.21.Hb,73.23.Ad 

I. INTRODUCTION 

The study of interacting fermions in one spatial dimension has attracted renewed attention in the recent past 
with the experimental realization of nano-structures and simulation of quasi-one dimensional environments 1] . Many 
theoretical works that seek to provide a basis for testing fundamental physics, as well as those that propose new 
scenarios for practical application of these ideas have emerg ed0. All these are based on the Luttinger liquid paradigm 
of HaldaneQ where homogeneous one dimensional systems are considered to be a specific kind of anomalous metal 
characterized by power law singularities in the Green function. However, this paradigm is valid mainly for short-range 
interactions. For a specific kind of long-range interaction (v(q) = 2e 2 /q 2 ), one of the authors have shown in an earlier 
workM that the system is characterized by essential singularities rather than power-law singularities. Also, in another 
work[5j, the formal theoretical framework upon which the Tomonaga-Luttinger picture is based, was critiqued. In 
particular, this approach fails to provide a proper description of the ground state of a one dimensional homogeneous 
Fermi system interacting via the gauge potential as has been pointed out Q . That the presence of Klein factors in 
this approach, renders the analysis of the three-wire junction problematic has been noted by Chamon et al @ and by 
Lai et al@- 

In this paper we show that fermions in one and higher dimensions may be conveniently bosonized by exponentiating 
the slow part of the field using current algebra. In particular, it is shown that no Klein factors need be invoked 1 . The 
right movers and left movers are shown to obey fermion commutation rules nonetheless. In one dimension, this fact 
is already known to the experts 8] 9j, however in higher dimensions it is probably less well known. Nevertheless, we 
present the one dimensional formalism even if there is some overlap with established methods, if only to motivate the 
generalization to higher dimensions. We also express the annihilation operator in momentum space directly in terms 
of a certain combination of number conserving products of Fermi fields known as sea-bosons |a| thereby completing 
the formalism introduced earlier 0- It is worth emphasizing that it is the slow part of the field operator that is easily 
expressed in terms of the bosons rather than the full field as was implied in one of the authors' earlier work 5j. The 
present approach, together with the action mentioned in the abstract, can be argued to be quite adequate in studying 
gauge interactions in higher dimensions as well. 



II. FORMALISM 



For a pedagogical introduction to bosonization in one dimension, there are many reviews available[9j- For the 
sake of brevity, we shall only discuss the important ideas. First we consider fermions in one spatial dimension. 



1 More precisely, just one global Klein factor suffices in all dimensions rather than one for each kind of mover|^], 0. 
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Define the right movers and left movers as usual (A <C kp is the momentum cutoff around the Fermi points), 
= -^zJ2k elkXc k F +kO{A- \k\) and -0 L (x) = J2 k e %kx c^ kF+k 6{k- \k\). These are slowly varying in space. In 
this theory we assume that in the end kp ^> A — > oo. Now we define the density fluctuation operator and its conjugate 
as follows ( with p and q unrestricted ), p q = J2 P c l+ q /2 c p-q/2 and X q = J2 P (^w) c P - q /2 c p+q/^- Tnese operators 
obey the following commutation relations ,[X q ,X q '] — [p q , p >] — and [X q ,p >] = i S q Since only terms with 

small q contribute ( the ip 's are slowly varying ) we have,[^(a;), p q ] — e lqx 4> r(x) ,[4> R,( x ) > X q ] = e~ lqx {jjj^ ^r(x) 
and [iJjl(x), p q ] = e tqx ipL(x), \>Pl(x), X q ] — — e~ lqx {jjj^j ''I'l(x). These may be exponentiated as follows. 

fef^e-^.^'e-^' 5 ?''^ (1) 



$ L (x) ^e-^^^e^^-^f" (2) 

The multiplicative factor, namely comes from the observation that, (^ip R (x)ipR(x)\ = (%l>^ L {x)ipL{x)^j = 5^. 

The Eq.(© and Eq. obey fermion commutation rules without the need for Klein factors for x 7^ x which is 
easy to verify. In the usual approach to bosonization, these formulae are regularized by introducing a decaying 
exponential 9] . The regularization is needed to reproduce the delta function in the fermion commutation rules, 
{ip R / L (x),ip R ^ L (x )} = S(x — x ). Here too, we adopt a similar but not quite the same approach. In this regard we 
postulate, 

^(*) = e-'E.^«e-E,^ ( * + 1 £ e -^_ pq + ] (3) 



That is, we add the density fluctuations of right/left movers to the (nearly infinite) mean density. This means that 
the commutation relations with p q and X q are now going to be invalid but the corrections are inverse in the number 
of right and left movers. Since A — > 00 we may write, 

^ fl ( x ) =e -*E,^e 4, "*, e -*E,*, T §=^/», e ^£^ e- <<! *e- |al/A [^+i^-«] (AY ( 5 ) 



where e — > + . It can be shown that Eq.© and Eq.© obey the following commutation relations, 

{iP R (x), i/> L (x)} = {^ R {x),^ L {x')} = {^ R {x),^ R {x')} = = ( 7 ) 



{^ R {x)^ R (x)} = ty L (x), Vl(x')} = S(x - x) (8) 

These nuances are not very important for the practical computations. In fact, we may use the unregularized fields for 
computing the propagators with finite spatial and temporal separations with impunity. We remarked that there are 
no Klein factors in the formalism. This refers to the need to invoke two different Klein factors, one for each kind of 
mover, in the usual approach 9J. In the present approach, there has to be an overall global Klein factor namely X q= Q 
that was defined in our earlier work0. We may now compute the propagators as follows. It can be shown that for 
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the free Fermi theory, (^ R (x ,t )ipR{x,t)j ~ ^ ^ ^ ^ py . One may also similarly study the interacting system 

with equal ease and obtain the right exponents for the Luttinger model. 

Next, we wish to generalize these ideas to higher dimensions. This is likely to be complicated since there is 
nothing called a right mover or left mover in two and more dimensions because the Fermi surface is no longer 
made of two points but rather, an uncountable infinity of them. The usual approach of Luther and HaldaneQ 
breaks up the Fermi surface into patches where the separation between the patches Akp ^> A ~ Aq. Not only 
is this contrived, it involves the introduction of one Klein factor in every radial direction. Here we show that 
there is no need for Klein factors at all (apart from the global one 7J). To this end we have to make the following 
analogy. Consider the following operators in one dimension, ipsiow{x) = -j= J2k eA k ^ kFSgn ^ x ct 9(A — \\k\ — fejr|) and 

ipdiff{%) = ttj; Sfc e t V°~ kFa9n V°)> x Cksgn(k) 9(A — \ \k\ — kp\). These are linear combinations of right and left movers 
which can be obtained from, 

ip(x,y) = ^^2e l ^- kFS9n ^ x e tS9n{k)v c k d(A - \\k\ - k F \) (9) 

as ip s iow{x) = ip{x, 0) ; ^diff{x) = —iip y (x, 0). By analogy we may generalize these ideas to more than one dimension 
as follows. We make the following prescription: sgn(k) — > k. Then, 

V>(x,y) = JL £y(k-fc,*).* e ife.y Ck 0( A -||k|-M) (10) 
v ^ k 

First, this object is a slowly varying function of x. In one dimension, a drastic simplification is possible which is 
not available in more than one dimension, namely, e isgn ^ v = 6{k)e iy + 9{-k)e~ i '« . In more than one dimension, all 
possible directions k are involved. It is this ^(x, y) which is a slowly varying function of x that is easily expressed in 
terms of the currents and densities rather than the full ip(x). This field ip(pc,y) is the higher dimensional analog of 
the one dimensional version where only two directions are involved and would correspond to linear combinations of 
right and left movers. Since this is a slowly varying function of x we can expect this to involve only pq and Xq for 
small q as is the case in one dimension. Define p q and Xq as in Ref. 5]. Then we have, 

^(x,y), Pq ] = -^£ e ^-^)-x e *y Ck _ q 0(AH|k|-M)«e~ i(Vy " q)(V,r ' x) ^(x,y) (11) 



[^(x,y),X q ] = _^£ e i( k -^)-* e *y 4,p- q / 2 Cp +q / 2 9(A- ||k| - k F \) « e^'*^'^ ^(x,y) 

These two rules may be exponentiated as follows. 

^>(x,y) =e-*^ e y y x « e E « e l^^J^ i?( X) y ) (13) 

Here F(x, y) is an integration constant independent of both Xq and pq. It obeys the constraint that VyF(x, y) = 

—F(x, y) since the unit vector k multiplies y. In practical computations, one computes the propagator of the 
interacting theory using the prescription in Eq. 1|13[) and multiplies and divides by the propagator of the free theory 
and in the division uses the one obtained from Ea. I|13fl and in the numerator uses the one obtained from elementary 
considerations. This eliminates F(x, y) and a closed expression for the full dynamical (albeit slow part) of the Green 
function may be written down in terms of current-current, current-density, and density-density correlations. The full 
propagator may be related to the current and density correlations as follows. 

(^(x\yVMx,y,i)) = e' E «(^) V ^^^ 
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e^^"^"' > «^^^<^-(*'W*)> ^t (x ' >y ' jt ' Wx>y>t)) (14) 



where -C •■■ ...> — < ... >o and the subscript < ... >o refers to the quantities for the noninteracting system. 

Next we wish to relate the fundamental field ^( x ) to the slow field if)(x, y). This exercise is essential in order to ensure 
a logical progression of ideas, also, it is the Green function of the fundamental field that is of physical significance. 
To accomplish this we may invert Ea JlOj l as follows. 

I f ^-tfor^ j = 1 6{k-\w\) (15) 

Since the vector y is dimensionless, we have to have an arbitrary (formally infinite ~ (kpL) d ) dimensionless 'volume' 
vq to ensure that the transforms with respect to y are properly inverted. Thus the fundamental field is, 

4>{X) = W ^ e *P-* e -i(|p|-**M>.* f QLe-V-Y^y) = - fcbcf dy E( X ,y;X) ^(x,y) (16) 

p J V Jvo v «0 J J VB 

where E(x, y; X) = y ^ p e »P- x e - l (p- fc F)p- x e - l P-y, In one and three spatial dimensions we have respectively, 



/ (fc F X-y).(X-x) \ 

k% cos { [x^x] ) 



E(x, y; A) « 2 cos{k F X - y) 5(X - x) ; £(x, y; X) » ^— J ^ (17) 

2n |X — x| 

In order for the one dimensional version of EcL ljlfiJI to be consistent with Eq.QJ and Eq.(J2J we have to ensure that 
F(x, y) = \J^- cos(y). One may now write down the full propagator of the fundamental field in terms of current and 
density correlations by replacing the right hand side of the equation below by the expression for the slow propagator 
in Ea. ljTH as follows. 

<^(X,t)^(X,t)>=\ fdx' I dy £(x',y';X') f dx f dy £(x, y; X) < ^(x, y, t)V(x' , y , t) > (18) 

We shall not dwell on this any further, for example, as a nontrivial application, we could calculate the full dynamical 
propagator of a two dimensional electron gas with long range interactions ( V(r) = Log(r) ). This is bound to 
lead to the conclusion that the system is a Luttinger liquid with a characteristic exponent, a result that will be 
more easily and elegantly shown using sea-bosons later on in this article. The current and density correlations may be 
conveniently calculated [at least at the RPA (random phase approximation) level for translationally invariant systems] 
using the Lagrangian approach that is conducive to the introduction of gauge fields [Tol|. To this end we derive an 
action functional for the free electron system in terms of the hydrodynamic variables in any number of dimensions 
at the RPA level using the sea-boson formalism. Then we argue using general methods, that the functional form of 
the general action should be similar to the RPA result. First, the derivation of the RPA form of the action in terms 
of hydrodynamic variables. We borrow the definition of the hydrodynamic variables in terms of sea-bosons from our 
early work0. 

x « = ™r- E — [^(q) + 4(-q)] ; = E^(-q) + 4(q)] ( 19 ) 

q k k 

We invert these using smearing functions. Let us postulate, 

A k (q) = -i r x (k, q)A q + T p (k, q)p_ q (20) 
The above is meant to be valid in an average sense and is not an operator identity. If we set, 

T x (k,q) = L4 k (q),4(q)] ; r p (k,q) = -J— ^ [v4k ( q)) 4( q)] (2 1) 

2Aoe q m 

then substituting Eq. i|2(J|) into Eq. (|19fl leads to an identity. We know from the sea-boson theory that the hamiltonian 
of the free Fermi theory is H = X)kq ^^(<l)^k(<l)' Thus the action in terms of sea-bosons 



is. 



s = [ l ^dtJ2 4(q) ( id t - ^) M<C ( 22 ) 



kq 
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Using Eq.JHI) in Eq.JHJ) we find 



s = [ *V[2r p (k,q)rx(k,q)] Pq d t A q - f d^r p (k, q ) 2 (^] PqP _ q - f d^r x (k,q) 2 A q A_ 

•/O ,,,, Jo ,,„ V "' I JO V "' I 

(23) 



k q " u k q " u k q 

This may be simplified to yield, 

-ifi r -i0 f-ifi 



S= f 1 dtV/, q fi,X q - f 1 dtY^wW-v- I ' dtJ^Noe^X^X-* (24) 
Jo Jo - /v o Jo 



q q 

2 



where eo = Y]^ 4Ar ^ e ^f(k) (^r) • I n one dimension eo = It appears that we have to use the one dimensional 
version of the energy scale eo in higher dimensions as well to reproduce the right static structure factor. The reason for 
this is similar to the argument that allows Schotte and Schotte .llj to use the Tomonga-Luttinger theory to analyze 
X-ray absorption in bulk metals in three dimensions (only the s-wave contributes). We now argue that this action is 
the RPA limit of a more general formula for the action in terms of hydrodynamic variables that will be postulated 
rather than derived. The claim is that the general action for the free Fermi theory is given in terms of the current 
and density variables as follows. 



P (vn) 2 



(VP) 2 
4p 



pd t u - v F ([p];x) 2m * p 1 i '->r,-. 

Here p and II are conjugate variables and J(x, t) — — p(x, t) VIl(x, t) is the current. Also p(x, t) = 
V S q n e_lq Xe_z "Vqn an d n(x, t) = Ylqn e* q ' x e s: ' l *Xq n where z n — 2im/(3 is the bosonic Matsubara frequency and 

J = f~ il3 dtjd d x. Here V F ([p];x) is a suitable functional of the density chosen so as to reproduce the properties of 
the free Fermi theory. To motivate the introduction of this new formalism we point out that in an earlier workp^ one 
of the authors had argued that the fermionic field (operator) may be expressed in terms of the conjugate variable as 
follows. 

V>(x, t) = e »A([p]; X ,t) e -.n(x,()y^ pjg) 

This means that the current operator may be written as follows J(x, t) = -p(x, ()VII(x, () + C([/)];x, t), where C 
depends on the A above. Such a general correspondence has already been suggested in the pioneering investigations 
of Rajagopal and Grest 12J. An examination of the commutator [Jj(x, t), Jj(x , t)} in terms of the Fermi fields shows 
that it is expressible as linear combinations of terms involving the currents themselves. By imposing this requirement 
on the bosonic representation for the current, we find that C = 0. This in turn means that A should be independent 
of x in order for current algebra to be respected. Unfortunately this conflicts with the requirement that A obey a 
certain recursion relation that was introduced in an earlier work 5] . Thus we take the point of view that A depends 
weakly on x so that the current is still given by J = — p\7H. However, the crucial A makes its appearance in the 
action as has been depicted in Ea. H25l) . The viewpoint advocated in this article is that we shall not be too insistent 
on verifying fermion commutation rules, rather this issue may be cleverly circumvented by forcing agreement with 
the correlation functions of the free theory and also by ensuring that the commutation rules of Fermi bilinears which 
are simpler, are recovered properly. Since Vf is uniquely determined by this approach, we may work backwards and 
determine A and then see if fermion commutation rules are indeed being obeyed. This is likely to be a formidable 
task and is therefore omitted. There are some subtleties associated with making the transition from the Hamiltonian 
description to the Lagrangian description. The principal of the these is the question of normal ordering. We note 
that the action in Ea.Q22|) is properly normal ordered. The hydrodynamic variables in Eq.JSBJ) have to be expanded 
in terms of the sea-bosons and terms have to be rearranged so that all the annihilation operators are to the extreme 
right. Then one de-quantizes Ea. (|26|l by demoting the various operators to c-number functions before using it to 
derive the action. The claim is that the resulting ordinary complex number valued function 'simulates' a Grassmann- 
valued field. The anticommuting property of the Grassmann variable is captured at the level of the propagator by 
the KMS(Kubo-Martin-Schwinger)-like boundary condition obeyed by the global Klein factor introduced in an earlier 
work|7j- As it is clear that normal ordering merely redefines A, the actual nature of the function A is determined by 
making contact with the free Fermi theory. To this end, we may expand the action in Ea. (|25|l in powers of density 
fluctuations. Retaining only the harmonic terms amounts to using the RPA. At this level we may write, 

A([p];x,t)« e^ x e^ 4 A qn p_ q ,_„ (27) 

q7 t0,n^0 
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where A q „ = - e ^ fl) and N° is the total number of electrons and tp is the Fermi energy. Therefore x) = pdt A. It 

can be shown that the action in Eq. (|25|l reproduces the right current and density correlations of the free Fermi theory 
at the RPA level provided we choose A to be given by Eq.JSZJ). In particular the choice in Ea. (|27|l substituted into 
Ea. (|25[) yields the RPA level action Ea. (|24ll derived earlier rigorously using sea-bosons. One may then contemplate 
adding more terms to this action such as density-density interactions, density-phonon interactions or current density- 
vector potential interactions to the action and expect to obtain nontrivial results for the one-particle properties. We 
may also make Eq. l|25[l plausible by pointing out that it implies a connection between the generating function of 
density correlations and current correlations. This is relegated to the appendix. 

In fact, we may go further and derive a general expression for the generating function of density-density correlations 
of interacting Fermi systems in terms of the corresponding correlation functions of the free Fermi theory. This exercise 
is straightforward to perform mainly because the action functional in Ea. (|25|) has a simple dependence on the conjugate 
variable II which may be integrated out. To this end, we may define the generating function for the density-density 
correlation functions as, 

Z([U]) = f D[p] f D[U] e lS e ^^ x U{x,t) P {*,t) ^ 



Correlation functions involving only the density variable may be obtained as appropriate function derivatives with 
respect to the source U, for example, < T [p(x,t)p(x ,t )] >= S 2 Z([U])/SU(x, t)SU(x , i )\u=o- The equation above ( 
Eci. (|28[) ) may be formally inverted and an expression for the action with the conjugate variable integrated out, may 
be written down in terms of the generating functional, 

D[IL] e lS = f D[U] e^fo'" 3 ' f d "* u' ( 2 9) 



One may now use the above expression to relate the generating function of the free theory with that of the interacting 
theory. Thus, 

Z ([U]) = ( d [p] [ D m e lSfr "e~i -fo' 113 dt S ddx J dd ' x ' ^(*-x')p(x>fM*'>*) e jo" ,/3 1 d " x u(x,t) P (x,t) 



= j D[U] J D[p] e~4^o W dtJ d d xj d d x V(x-x')p(x,t)p(x',t) e J ^ I/3 / d d x (t/(x,i)-C/'(x,t))p(x,t)^ ee ^jj'j^ 

We may now decompose the density variable into the various Fourier modes and perform the integration over the 
density variable and write, 

Z([U]) = J D[U] e E ^°" mkV^-K.nW-v.-n-uL^J Zfree (p'^ (31 ) 

where Vq = J d d x V(x) e~ tq x and Uq n — y J Q ^ dt J d d x J7(x, i)e~ jq - x e~ 2,,t . The generating functional of the free 
theory may be written as follows. 

Z free {[U}) = e* s ^" <p^p-«,-~>o u^U-^-n e AE WiB / <^»/» Q ',„'P-,-,',-„-„'>o t^tV in 'tf_,_,',-„-„'... (32) 

Thus we may rightly suspect that retaining only the gaussian terms in the above generating functional gives us 
the density-density correlation in the sense of RPA. Corrections to RPA are obtained by retaining the higher order 
terms involving three-body density correlations. Unfortunately, the current-current correlations are not so simple. 
This is because the current is a nonlinear combination of the density and conjugate variables. However, one may 
contemplate reverting to the Hamiltonian description and invoking the equation of continuity to relate the current- 
current correlations with the density-density correlations. 

Next we wish to show how to express the annihilation operator in momentum space c p directly in terms of sea- 
bosons[3 in contrast to the exercise just completed, where we expressed only the slow part of the field in terms of 
the current and densities. To be sure, the exercise in this section is also valid only in the RPA-sense and hence only 
provides the asymptotics. Define n q (k) = c k_ q/ / 2 c k+q/2 f° r Q 7^ 0. Hence we may write, 

[C p , n q (k)] = <5 p k_ q /2 Cp+q = <5p,k-q/2 ^(p) C p (33) 

where T q (p) = exp(q ■ V p ) is the translation operator. Also, 

[n q (k),n_ q '(k')] = [4- q / 2 c k+q/2, c I' +q ' /2 c k'-q'/2] ~ ("f (k - q/2) - n F (k + q/2) ) Vk'^q.q' =sgn(k.q) <5 k , k <<5 q , q < 

(34) 
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The last two approximations are equivalent to the random phase approximation. The commutation rule in Ea. H33|) 
may be exponentiated as follows, 



- p-E q s9"(p-q+q 2 /2) «_ q (p+q/2) T q (p) 



c p = e 



/(p) (35) 



where /(p) is an integration constant independent of n q (k), the exact nature of which will be discussed subsequently. 
To be sure there are additional terms caused by the fact that the translational operator also acts on these commutators 
and so on. The reasons why these issues are not important are relegated to future publications. Suffice it to say that 
this is consistent with the RPA. Now we wish to evaluate /(p). To this end we make the following surmise0, 



/(pl^-^WM) (36) 

Here no(p) = c P c p an d fo(p) is canonically conjugate to n (p), in other words, [-Po(p), ^o(p )] = iS p >. This ensures 
that [c p ,rt (p )] = c p <5 p , P '- We have also set [n q (k), n (p)] = 0. In practical computations, we are going to set 
rio(p) ~ np(p) = 6(kp — |p|). This means we have to ensure that -Po(p) is formally infinite in order for it to be a 
conjugate to no(p). Thus we redefine, 

/(p) = e- iN °^ n F (p) (37) 
Here N a is a large quantity (total number of particles) and £ is an arbitrary c-number function. This means it obeys, 



Lt^e-^^e^'^S^ (38) 

Unfortunately this prescription leaves out an essential part of the dynamics. For example we expect that for the free 
theory, c p (i) = e~ 4£p *c p (0). This is not obeyed unless we retain Po(p) as an operator conjugate to ^o(p) 2 ■ Thus we will 
have to put the rapidly varying exponential in by hand during practical computations. More precisely, we are going to 
multiply and divide by the free propagator and use the bosonized version in the denominator and the one obtained from 
elementary considerations in the numerator thereby rendering this issue moot. In terms of bosonic operators defined in 
our earlier work @, n q (k) = A k (q) + 4k(-q). Therefore, sgn(p.q + q 2 /2)n_ q (p + q/2) = -A p+q / 2 (-q) + A p+q/2 (q). 
In other words finally we obtain, 

Cp = e£«»(W(-i)-^(i)) ^(p) e -w°«P) nF(p) (39) 

By expanding the exponential, it can be easily verified that this formula is consistent with the definition n q (k) = 
c k:-q/2 Ck +q/2- Notice that there are no Klein factors in this formula either. Practical applications of the above 
formula will be relegated to future publications. One may contemplate using Ea. H39|l to compute the one-particle 
Green function of any Fermi system. The results are valid at the level of RPA. In fact in an earlier workQ, one of 
the authors has highlighted the novel non-Fermi liquid character of a system described by long-range interactions in 
one dimension. In the following section, we study the same system in two spatial dimensions rather than in one. Here 
too we find that the system is an anomalous metal ( or an insulator ) but of the usual power-law Luttinger kind. 



III. METAL-INSULATOR QUANTUM PHASE TRANSITION IN TWO DIMENSIONS 



Now we wish to turn our attention to the properties of a system of spinless electrons mutually interacting via a 
gauge potential i> q = 2-7re 2 /g 2 in two spatial dimensions. In real space this corresponds to a Logir) potential. Theories 
of novel non-Fermi liquid properties in higher dimensions are rare one such study for example is by Gori-Giorgi and 
Ziesche^3|. We consider the following hamiltonian in the sea-boson language 0. 

H = E 0£) 4(qMk(<0 + E w Et^(-q) + 4(q)P k ' (q) + 4 (-<0] (*o) 

kq ^ m ' q^O kk' 



2 Since n q (k) varies slowly in time we may ignore it for this aspect of the discussion. 
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This represents electrons mutually interacting via a two body potential u q such that only small momentum transfer 
terms are important (forward scattering only) and k and k are close to the Fermi surface. Such a constraint is realized 
by the choice v q = 2ne 2 /q 2 since the potential is singular for small q and negligible for larger q. This hamiltonian 
may be easily diagonalized via a Bogoliubov transformation. We may borrow the results of Ando et al[l4j for the 
dielectric function of a two dimensional electron gas to arrive at the following formula for the dispersion of collective 

modes : w c (?) ~ "7Sr= + I — 3fcF 5 ) q 2 ■ The smallness of q is governed by a cutoff A which is the value of q for which 

V2-TO y4\/2em2 J 

the second term is of the same order as the first. This means A w -^e^/m. The momentum distribution of this 
system at zero temperature may be evaluated using the sea-boson formulas^. 

< n k >= 1[1 + e - 25 -( k >]n F (k) + i[l - e - 2S -( k '](l - n F (k)) (41) 

where S% (k) = £ q < 4- q/2 (<l)^k- q /2(q) > and S% (k) = £ q < 4 +q / 2 (<l) A k+q/2 (q) > and n F (k) = #(fc F - |k|). 
In the q summation we retain only the most singular parts and integrate only upto the cutoff A. We omit the rest of 
the details regarding computation of the boson occupations. When this is done we obtain the following formula for 
the momentum distribution. 

fcJMiy sgn(k F -M) (42) 

where 7 = -p- —£= — ^rwr- and A = e~ * A and uin — ^7= . Thus the system an insulator of the usual Luttinger 
type for ujq > 2Ef and an anomalous metal with a residual Fermi surface for luq < 2Ef- We may now compute the 
dynamical Green function using the methods outlined in the earlier sections, but this is not going to yield any new 
information, since the important attributes namely the estimate of the plasmon energy and the anomalous exponent 
have already been highlighted. It would also be interesting to see if the usual Coulomb interaction in three spatial 
dimensions also exhibits some kind of metal insulator transition at absolute zero temperature as a function of density. 
Unfortunately in order to study this reliably we have to go beyond the RPA and this calculation is beyond the scope 
of this preliminary work that focuses mainly on formalism. 



< rik >= 



IV. CONCLUSIONS 



To conclude, we have shown that it is possible to bosonize fermions in arbitrary dimensions and relate bosonization 
with quantum hydrodynamics. We have been able to extract the dynamical one-particle Green function in terms of 
the current and density correlation functions valid at the level of RPA. We have expressed the fermion annihilation 
operator in momentum space directly in terms of bosonic operators called sea-bosons. We have also derived a general 
formula for the generating function of density correlations of systems interacting via two-body forces in terms of the 
corresponding quantity for the free theory. Finally we have evaluated the momentum distribution of a two-dimensional 
anomalous metal ( or insulator ) with the electrons interacting via a long range potential and extracted the anomalous 
exponent. 



V. APPENDIX 



The connection between the current-current correlations and density-density correlation implied by the equation of 
continuity is recovered by the action Ea. (|25|l postulated in the main text. This may be verified at the level of RPA 
in any number of dimensions as follows. Define, 

/ P (vn)^+< v /> 2 , 
if pd t Il~pd t A 2m 4 " +/ Up 



Z d [U] = / D[p]D[U] e V / (43) 



Z C [V] = / D[p}D[H] 



i f pdtYl-pdtA-- 



-f J-VV 



(44) 
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Then < T[p(x,i)p(x',*')] >= , y(x ffiff xV) and < T[(V X • J(x, i))(V x , • J(x', i'))] >= t >y We may 

relate Z c and Zd as follows. 

Z C [V] = Z d [~ (W) 2 + m d t V] (45) 
The validity of this may be verified at the level of RPA. At this level we may write, 

Z d [u] = e !/< T M x >*)p(xV)]> U(x,t)U(x.',t) ^ 

This implies, 

<T[V X -J(x,t) V x /-J(x',t')]>= J <T[p(y,T)p(y',T')}> (m 8(y - x)d T 6(T - t))(m S(y' - x)d T ,5{T - t)) 

= m 2 <T[d t p(x,t)d t ,p(x',t)}> (47) 
The validity of the above equation follows from the equation of continuity. 

VI. DIRECTOR'S CUT 

In this section, we provide a new proof of the irrotational nature of the velocity operator for fermions. This makes 
the above paper logically complete. Unfortunately this proof is not there in the published version. 

Lemma: If Vp x J — pV x J = then there exists a II such that J = — pVII. 

Proof: First we note that [J(x), p(x)] = 0, hence we may write j without ambiguity. 

Vx- = ivxJ - 4rVp x J = 4r[pV x J-Vpx J1 = 
P P P 2 P 2 

Hence j = — VIT for some II. Hence J = — p\7U. Now wc prove that Vp x J — pV x J = for (free) fermions. We do 
this using the first quantized definitions of current and density operator and using real space fermion wavefunctions. 

N 
3=1 

. N N 

3=1 3=1 
N . N N N 

Vp(r) x J(r) = r 3 )} x £[V<5(f - r p )} £[V<J(r - fj)] x i£ 6(f r p )V p 

3=1 P=l 3 = 1 P=l 

N N 

p(r)V x J(r) = -iYj^r-r^YsNSir-rp)] x V p 

3=1 P=l 

N.N N 

V(r) - Yl [W~ fj)] x r p )] - i ]T 5 if~ r P )[VS(f - r})] x V p + i ]T <5(f - ?j)[VS(r- f p )] x V p 

J,P=1 3,P=1 3,P=1 
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. N 



V W = o E rj)mr- r P )\ - 8(r- f p )[VS(f- **•)]) * (V p - V,) 



Let, C(rj,r p ) = (S(r — fj)[V5(f — r p )\ — 5(r — f p )[V5(f — fj-)]) x (V p — Vj) and f(fj,r p ) be some function : 



y dr p C{rj,r p )f{rj,r p ) = J dr j df p {S(f-r j )[VS(f-r p )]-S{r-f p )[V5(f-f j )])x(f 2 (r j ,^ 



)) 



= ^ dr p [V<5(f - r p )] x /b(f,r p ) - J dfj[V6(f- fj)] x ^(rj,r) 



^ <F P [V<5(r-r p )] x h{r,r p )-J dr 3 [VS(r - r})] x /^r) 



= e/JK (h,J;2,K(r,r) - fi,j-,2,K(r,rf) - J^e/ e/ JK (^f 2 ,j-i,K(r, r) - f\,j-i,K{r,r)} 



fl,J;2,K(r,r) 



_d d_ 

dxj dy K 



/(X,Y) 



X=Y=f 



_9 9_ 

dxj dy K 



/(Y,X) 



sf2,J;l,K(f,f) 



X=Y=f 



h,j;i,K{r,r) 



d d 



dyj dy K 



/(X,Y) 



= s 



X=Y=r 



_9 a_ 



/(Y,X) 



= s/i,j ; i,Ar(r,r) 



X=Y=r 



••• = (f2,j;2,K(r,r) - fi,j-2,K(r,r)) - {h,j-i, K {r,r) - /i,j ; i,j<-(r,r)) 

= (s + - (s + l)/2,j;i,A-(r,r) 

For s = — 1 then these terms are zero. 
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